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where the sum is performed over non-strange hadronic fi- 
nal states hy^A with angular momentum J = 1. Ny/A 
is the number of taus that decay into the hadron hy/A 
plus neutrino, Wy/A denotes the appropriate weight of 
the hadronic mode to the vector or axial-vector current, 
'5'ew = 1-0194 is an electroweak correction term [||] and 
iKdP = 0.9477 ± 0.0016 is the squared CKM weak mix- 



ing matrix element | 

Within the framework of QCD weighted integrals over 
the spectral functions (so called moments) have been cal- 
culated lUT 
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Abstract 

It has been shown in recent analyses by ALEPH |jl|] and 
OPAL [|[] that precision QCD tests are possible with 
hadronic t decays by comparing spectral moments of the 
hadronic decay ratio of the t with QCD calculations. In 
principle e+e^ data can be used in a similar manner by 
evaluating spectral moments of R. The current e+e" data 
is compared with the OPAL t data and a prediction is made 
on the achievable accuracy of QCD tests with the projected 
precision ofPEP-N 

1 INTRODUCTION 

The T lepton is the only lepton heavy enough to decay into 
hadrons. The observed spectra of the non-strange hadrons 
in s, where y^ is the mass of the final state hadronic sys- 
tem, give the non-strange spectral functions. Decays with 
an even number of pions in the final state belong to the 
vector current while decays with an odd number of pi- 
ons belong to the axial-vector current. A comparison of 
weighted integrals over the spectral functions of the vec- 
tor and axial- vector current with QCD predictions can give 
fundamental parameters of the theory [Q-H], including the 
strong coupling constant a^. Different integrals (moments) 
are used to measure power corrections of non-perturbative 
origin and the strong coupling simultaneously, thus reduc- 
ing the theoretical uncertainties on as connected with the 
non-perturbative terms. 

The QCD tests can also be performed for energy scales 
smaller than nir, if the integrals over the spectral functions 
are performed with an upper integration limit rn^, — sq < 
and replacing iri^ by m^, in the integrals [y thus cre- 
ating the spectral moments of a hypothetical r' lepton with 
a mass below rrir- 

In the same manner e"'"e~ annihilation into hadrons can 
be used to extract the spectral function of the vector-current 
and compare with theory by means of its spectral moments. 

2 THEORY 

2. 1 Hadronic r decays 

The vector v{s) and axial-vector a{s) spectral functions 
(the absorptive parts of the vacuum polarization correla- 
tors) are given by the spectra in s of the final state hadrons, 
normalized to the branching ratios B and corrected for the 
phase-space; 



where the scalar spectral function v^{s) vanishes, since no 
scalar particle has been observed in r decays, while the 
pseudo-scalar spectral function a"(s) is given by the pion 
pole, assuming that the pion is the only pseudo-scalar final- 
state in non-strange r decays: 
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The moments are used to compare the experiment 
with theory. In what follows, ten moments for kl = 
00, 10, 11, 12, 13 for V and A ai-e used. The first mo- 
ments R'^'^y/A '■^^ normalized decay rates of the 
T into vector and axial-vector mesons. In the naive parton 
model these two rates are identical and add up to the num- 
ber of colors. Since only non-strange currents are consid- 
ered in this work the naive expectation has to be multiplied 
by I Kid p. Including the perturbative and non-perturbative 
contributions, equation (0) is usually written as [pi]]: 
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where 6^^^^^ are perturbative QCD corrections (w 20 % for 
kl = 00) and the (5yy^ terms are the so-called power cor- 
rections (« 1 % for kl = 00). 



The perturbative term Jp^j.^ is known to third order in 
as and partly known to fourth order in a^. For kl = 00 
the Contour Improved Perturbation Theory (CIPT) resuh 
is Pi: 
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where the Ki = 1,^2 = 1.63982, /fa = 6.37101 ai-e 
known Jl^-p^ and = 25 ± 50 has been estimated Jill , 
|l7|,|l8|]. The Taylor-expansion in as{m'^) of the CIPT re- 
sult is [|7l|n|]: 



"pert 



-5.2023a^ + 26.366a^ + (78.003+X4)as, (6) 



with as = as(m^)/7r. It is referred to as Fixed Order 
Perturbation Theory (FOPT). The third method considered 
in this paper resums the leading term of the /3-function to 
all orders in as by inserting so-called Renormalon Chains 
(RCPT) [[Ill-Ill. The fixed-order corrected version (up to 
the third order in as) quoted in the lower portion of table 6 
in reference [ ^0| ] is used in the fit. 

The power corrections (Jy/A '■^^ framework of the Op- 
erator Product Expansion (OPE) [ ^3| ] are proportional to 
. The dimension D = 2 terms are mass-corrections 
of perturbative origin [0,|ll[l and are small for non-strange 
r decays. Corrections of higher dimension are of non- 
perturbative origin, absorbing the long-distance dynamics 
into vacuum matrix elements [^^^-^]. In contrast to the 
perturbative part the power corrections differ for the vector 
and the axial-vector currents. 

If one neglects the small s-dependence of the power cor- 
rections, the <5y^^ terms can be expressed for all kl values 
by a product of the power correction for kl = 00 and a sim- 
ple integral over the fc^-dependent weight-functions [pT|]. 
This approach is used for the dimension D = 6 and D = 8 
terms, taking Sy^^" as free parameters. For the dimension 
D = 2 and = 4 terms the full s-dependence is taken into 
account for the theoretical description of the moments [pT|]. 
The least precisely known D — 4 parameter, the gluon con- 
densate which is known only to 50 % [Q] is also taken as a 
free parameter in the fit, while the D = 2 term is calculated 
from the quark masses and the strong coupling. The quark- 
masses and -condensates needed to complete the 13 = 2, 4 
terms are taken from [0] . 

2.2 e"'"e~ annihilation into hadrons 
The ratio R^+c- is defined as: 

Rc+c-{s) — ^^^2 '''c+c- ^hadrons i^) 

= UnlmWis) ^6v''{s), (7) 

with v'^{s) « X]/=u d s Q}'^{^) being the vector spectral 
function with (except for the isoscalar contributions and the 



charge dependent factor) similar properties as v{s) in ([p 
and (^ and Q / denoting the charge of the quark flavor /. 

In massless perturbative QCD J^f Q'f''j{s) and v''{s) are 
identical and most conveniently expressed in form of the 
Adler-function f^: 
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with (for three flavors) the same Kn as in (||). 

The mass corrections to the Adler-function are given 
by M: 
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where the running quark masses mf{s) are calculated from 
scale-invariant mass parameters and an evolution equation 
which is know to four loops [^9|]. 

Finally the non-perturbative part of the Adler function 
reads [0]: 
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with the dimension 4 contributions from the gluon conden- 
sate {as/nGG) and the quark condensates {mfipf^'f) and 
the dimension 6 and 8 operators (Og) and (Og). 

The total Adler function is simply the sum of the per- 
turbative part (^, the mass corrections (^ and the non- 
perturbative part (|Tc|): 

D{-s) ^ Dpi-s) + D,^,ss{-s) + D^p{-s). (11) 

Moments of similar form to (^) can be defined for e+e~ 
data M]: 



Ami 



In order to guarantee the validity of the OPE the end- 
point in the integrals should be suppressed by at least two 
powers of s and the moments chosen here are therefore 
kl = 20,30,31,32,33. As in r decays the moments can 



be rewritten as integrals around the circle |s| — in the 
complex s-plane: 



/iS„-(»o) = 6»i^ (13) 



|s|=so 



and further written as integrals over the Adler function after 
integrating by parts: 
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In the QCD fits the dimension D = Q and = 8 terms 
are fitted using (Oe.s) as free parameters. Note that this 
differs from the r data fits where the contributions to the 
moment kl — 00 are used as free parameters. For the di- 
mension D = 2 and D — 4 terms the full s-dependence 
is taken into account for the theoretical description of the 
moments as described above. In contrast to the non-strange 
T data where the least precisely known D = 4 parameter is 
the gluon condensate the largest uncertainty for the e+e^ 
data comes from the strange quark mass ms(l GeV) which 
is used as a free parameter for the D — 2 and D — 4 terms. 
The other quark-masses and -condensates needed to com- 
plete the I? = 2, 4 terms are taken from [Q]. 

3 RESULTS OF QCD FITS TO r DATA 

The fit results of QCD parameters to the ten moments 
K\/A for = 10' 11' 12' 13 as reported in |@,|3l|] 
are shown in Tab. |l|. Only the CIPT fits are presented here, 
since the focus is the stability of the QCD fits and the two 
other approaches lead to similar results within the theoreti- 
cal uncertainties. The various fits demonstrate the stability 
of the method and that the perturbative parameter can 
be measured together with the non-perturbative parts even 
if most of the non-perturbative parts do not cancel as it is 
the case in the first three fits. 

Note that the statistical error on as in the last two fits 
is different only because the last fit uses additional infor- 
mation from the r lifetime and the leptonic branching ra- 
tios [gl. 

The extension to energy scales s < ni^ is demonstrated 
in Fig. Here the moment v+Ai^o = "^^0 shown 
together with the fit results to the three different theories 
FOPT, CIPT and RCPT for the perturbative part plus non- 
perturbative parts. 

The error on the theoretical curves (or alternatively their 
spread) indicates that the OPE is applicable to values as 
low as 1.5 GcV^ which can be used to extract the running 
of the strong coupling constant from t data alone. 

The three considered theories differ in the treatment of 
unknown higher order terms which is partially taken into 
account as theoretical uncertainty by means of the error on 
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Figure 1: R^,'^^(so)fora hypothetical t' lepton versus the 
upper integration limit sq — rii^,. 

K4. As a conservative approach one could take the average 
of all three results and their spread as additional theoretical 
uncertainty. From the results in |^] one obtains therefore: 

as{ml) ^ 0.326 ± 0.007cxp ± 0.022thoo, (15) 

which translates to 

as(m|) = 0.1194 ± O.OOOScxp ± 0.0027thco (16) 

atmz 91.188 GeV. 

4 AVERAGING OF e+e" DATA 

Exclusive c+c^ channels up to a CMS energy of 
2.2 GeV [p2[-|60|] are combined in this analysis to calculate 
the total hadronic cross section. The narrow ui and (p reso- 
nances are excluded from the exclusive channels and added 
to the total cross section as relativistic Breit-Wigner curves 
with s-dependent widths [ p^ |6lt l . The weights of individ- 
ual channels and the treatment of systematic errors is taken 
from [ ^0| ] but the statistical errors and the final result are 
obtained in a slightly different manner: 

1 . a common equidistant binning in energy is chosen for 
all channels (N = 300 bins from £'i„i„ = GeV 
to -Emax = \/5GeV) and for every experiment and 
channel a histogram is filled with weighted averages 
of cross section measurements falling in the same bin: 



{j\NEj/E„-,^^e[i-l.i]} 
{j\NE,/E^^^G[x-l,x]} 



, i = l...N, (17) 



Table 1: Comparison of the QCDfit results ( CIPT) to the moments of vector (axial-vector) current, the simultaneous fit of 
all moments for both currents, and the moments for the sum of both currents /|37|/. The quoted errors contain statistical 
errors only. 
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where the squared errors Af are given by the inverse 
of the denominator in ( 17 ). 
2. gaps at fc ~ i ■ ■ ■ j between bins i and j of indi- 
vidual measurements are interpolated with the trape- 
zoidal rule: 



After this interpolation and averaging procedure the total 
hadronic cross section is given by the sum of the exclusive 
channels as in [|o[. The resulting values for R^+c- 3xq 
shown in Fig. 0. 



dk = Ckdi + {1 - Ck)dj , Cfc 
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3. the statistical errors are interpolated with the same 
procedure: 



Afe = CfeA, + (1 -Cfc)Aj-, 



(19) 



and are therefore larger compared to Gaussian propa- 
gated errors by the factor: 



rk = 



CfcAj + (1 - Cfc)Aj 
CfcA, e (1 - Ck)Aj 



(20) 



4. the correlation matrix for the interpolated parts is 
given by: 



V'stat 
kl 



TkTl {ckcAl + (1 - Cfe)(l - c^Af) . (21) 



5. the systematic errors are also interpolated according 
to the trapezoidal rule (|l9|) and all systematic errors 
are assumed to be 100 % correlated. 

6. the final error matrix is given by the sum of both: 
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7. the averaged cross section over all experiments in 
bin i for one exclusive channel is obtained from the 
weighted mean of all measurements and its error by 
usual error propagation. In case of inconsistent sets of 
data points for a particular bin its final error is scaled 

by the ratio S — \Jx^ /xls%' if S* > 1, where is 
the calculated for the individual experiments being 
consistent with the average and is the value 
corresponding to a confidence level of 68 %. 
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Figure 2: ^- {s) from exclusive channels. The dark gray 
band denotes statistical errors; the light gray band shows 
the sum of systematic and statistical errors. The dark solid 
line shows the relative uncertainty; the small dark error 
band beyond s = 2 GeV^ shows the projected uncertainty 
after 5 years of PEP-N j^. The narrow uj and (/) reso- 
nances and the massless QCD prediction are also indicated 
by solid lines. 



RESULTS OF QCD FITS TO e+e 
DATA 



i?^'i^_(4GeV^) for kl 
nven in Tab. ^ Their correla- 
I The result of the CIPT fit to these 
moments is given in Tab. The correlations of the fit 
parameters are shown in Tab. 



The five moments 
20,30,31,32,33 ai-e : 
tions are given in Tab. t 



Table 2: Moments i?* 



for 



kl 

20, 30, 31, 32, 33 /rom exclusive e"^e~ data. The er- 
rors include statistical and systematic uncertainties; 
projected errors after 5 years running of PEP-N and 
theoretical uncertainties are also given. 
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Table 3: Correlations of the moments R'^i^^ (4 GeV^)for 
kl — 20, 30, 31, 32, 33 in percent. The values correspond 
to the quadratic sum of experimental and theoretical er- 
rors. 
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Table 4: Results from the QCD fits to the moments 
R'^i^^iAGeV^) for kl = 20,30,31,32,33 f mm exclu- 
sive e~^e~ data. The errors include statistical and system- 
atic uncertainties; projected errors after 5 years running of 
PEP-N and theoretical uncertainties are also given. 
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The value of the strong coupling at 4 GeV corresponds 



to 



a,{ml) = 0.117 ± O.OOScxp ± 0.002th, 



(23) 



at niz = 91.188 GeV in good agreement with the t re- 
sult 16 but with a larger experimental uncertainty. Figure ^ 



Table 5: Correlations of the QCD fit parameters to the mo- 
ments (4 GeV^) for kl = 20, 30, 31, 32, 33 in per- 
cent. The values correspond to the quadratic sum of exper- 
imental and theoretical errors. 
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Figure 3: ihe moments n,^+^- (sq) versus the 

integration limit sq. The shaded bands show the experi- 
mental moments with statistical and systematic uncertain- 
ties from top to bottom in the order kl — 20, 30, 31, 32, 33. 
The dotted curves denote central values and ±lcr ranges 
for the theoretical expectations (CIPT) using the fit values 
to the moments at sq = A GeV^ as input. 



shows the moments of i?c+o- as a function of the up- 
per integration limit together with theoretical predictions 
using the fitted values at so = 4 GcV^ as input. Rea- 
sonable agreement between the extrapolated and the mea- 
sured moments is observed for all five moments down to 
So ~ 2.5 GeV^. The first three moments agree within er- 
rors with the extrapolated results down to sq ~ 1-5 GeV^. 

6 CONCLUSIONS 

The comparison of QCD fits to different currents (V, A, 
V and A, Vh-A) in non-strange hadronic r decays probes 
the stability of the OPE at low energy scales. The good 
agreement of perturbative and non-perturbative parameters 
among these fits demonstrates that QCD fits can reliably 
be performed at these low energy scales. The use of spec- 



tral moments is therefore extended to e+e~ hadrons 
data in order to fit QCD parameters at 4GcV'^. While 
the extraction of the strong coupling as is not competi- 
tive in terms of experimental uncertainties with the mea- 
surement from T data it provides an important cross check 
and probes an energy region not accessible with t de- 
cays. Furthermore the sensitivity to the mass of the strange 
quark can be used to constrain this parameter. A value of 
TOs(l GeV) = (220 ± 36 ± 59) MeV has been observed. 
The experimental uncertainty could be reduced by 30 % af- 
ter 5 years of running with PEP-N. 
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